Adiabatic preparation of many-body states in optical lattices 
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We analyze a technique for the preparation of low entropy many body states of atoms in optical 
lattices based on adiabatic passage. In particular, we show that this method allows preparation of 
strongly correlated states as stable highest energy states of Hamiltonians that have trivial ground 
states. As an example, we analyze the generation of antiferromagnetically ordered states by adiabatic 
change of a staggered field acting on the spins of bosonic atoms with ferromagnetic interactions. 
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Preparation and characterization of strongly correlated 
phases such as magnetically ordered states is one the 
more intriguing directions in the field of ultracold atoms 
[T]. In principle, such states can be prepared by start- 
ing with an ultra-cold atomic gas, turning on the opti- 
cal lattice, and reaching the magnetically ordered Mott 
state [2H1]- The direct preparation of such states is, how- 
ever, challenging as the energy scale of magnetic super- 
exchange interactions is quite small. An attractive al- 
ternative approach is to adiabatically prepare the de- 
sired states starting from a more easily produced ini- 
tial state [5l411|. In this paper, we demonstrate that 
adiabatic preparation of strongly correlated many-body 
states such as antiferromagnetically (AF) ordered states 
is indeed feasible under current experimental conditions. 
Furthermore, we point out that this allows for realiz- 
ing interesting quantum phases by adiabatic preparation 
of the highest energy states of underlying Hamiltonians, 
which may not have interesting ground states. This ap- 
proach utilizes a unique feature of systems of ultracold 
atoms, namely their nearly perfect isolation from the en- 
vironment. Such isolated systems with bounded energy 
spectra can have stable, low-entropy states near the max- 
imum total energy [T2II3]. The present work extends this 
concept to strongly interacting many-body systems. Un- 
like weakly interacting systems, the highest energy states 
of strongly correlated spin systems can be substantially 
different from and are often are more intriguing than the 
corresponding ground states. A dramatic example of long 
lived high energy metastable states was recently provided 
by the experimental demonstration of repulsively bound 
pairs [14]. 

Before proceeding we note that in Rcf. 7J a related 
approach to adiabatic preparation of gapless AF ordered 
states was investigated starting from a state with an ex- 
ternally constructed staggered magnetization. In this 
scheme, the Hamiltonian is gapless throughout the prepa- 
ration stage, and it is hard to maintain adiabaticty. In 
our approach the state is protected by a large gap dur- 
ing most of the evolution, which makes the preparation 
much more robust. Furthermore the possibility of study- 
ing the highest energy state was also noted in Ref. [5]. 



We provide general requirements for the applicability of 
this method based on time reversal symmetry. 

The main focus of this work is the Heisenberg Hamil- 
tonian for spin-1/2 particles on a lattice 

H = -J^S,S,. (1) 

This model with J > describes spin exchange interac- 
tions for two component Bose mixtures in the Mott state, 
when the scattering length is independent of the spin 
states [8]. This condition is realized for the experimen- 
tally relevant case of Kb atoms away from a Feshbach res- 
onance. The ground state of this model is ferromagnetic 
and not very interesting. On the other hand, the high- 
est energy state is AF with quantum fluctuations playing 
an important role. This can be understood by observing 
that the highest energy state of Ti. is the ground state of 
%' = Another interesting example was studied in 

Ref. [5 where it was shown that by using the highest 
energy state one can study the critical region as well as 
Haldane's phase for spin one systems in one dimension. 
We will also discuss how to realize a frustrated system as 
the highest energy state of spinless bosons. Frustrated 
systems have previously only been discussed for the spin 
state of ultra cold atoms [IS] . With spinless bosons one 
avoids the slow time scale associated with the spin dy- 
namics making it much less demanding experimentally. 

To be concrete we discuss the preparation of an AF 
state using two spin states of Rb atoms, focusing on an 
approach utilizing experimental tools similar to ones de- 
veloped recently at NIST [TTj- The NIST experiments 
combine magnetic trapping of two spin states of Rb, spin 
independent optical lattice, and the possibility to apply 
spin dependent staggered potential. 

The procedure is as follows: we start with a single 
component Bose gas in an optical lattice in the Mott 
state and apply a staggered magnetic field h{t). The 
effective Hamiltonian is then of the form 

n{t) = -jJ2s,s, + h{t) (y, St - J2 , (2) 
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with A and B being two sublattices. If h{t) ^ J mi- 
crowave or rf coupling can selectively spin-flip every other 
atom to a different spin state [TT] producing a two com- 
ponent Bose mixture with AF order, e.g, 
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We assume that the initial staggered field is large h{t) ^ 
J, so that the ground or highest excited states are unique 
and have large gaps to spin excitations, but that it is 
smaller than the onsite interaction h{t) < U so as to 
not destroy the Mott order. Preparation of the ground 
or highest energy states is determined by which set of 
spins are flipped, resulting in AF order either aligned or 
anti-aligned with the staggered field. 

For the anti-aligned configuration (|3| the system is in 
the highest energy state of the spin Hamiltonian ([2), 
which is gapped from even higher energy particle-hole 
excitations in the Mott sate. As h{t) is ramped down to 
zero adiabatically, the system should stay at the highest 
energy state of the instantaneous Hamiltonian. When 
h{t) = 0, we find a system described by the Hamiltonian 
([T]) near its highest energy state. 

Our method provides practical advantages for realiz- 
ing ground states of interesting spin models as compared 
with direct preparation by loading into an optical lat- 
tice. First it allows to study the ground states of "H, the 
natural Hamiltonian realized in the system, as well as 
those of —T-L. Second, the loading of the atoms into the 
lattice does not require adiabaticity with respect to the 
weak inherent spin exchange interaction J. The latter 
is only required for the transformation from a state ([3| 
which is relatively close to the desired final state, and 
this evolution can therefore be rather fast. 

As an example of the second point consider the natu- 
ral realization of AF states by Fermions in an optical lat- 
tice. One of the main challenges for direct preparation by 
loading into an optical lattice consists in maintaining adi- 
abaticity with respect to the effective spin Hamiltonian 
as the lattice potential is raised. For the present proce- 
dure, on the other hand, we only require that initially a 
spin-polarized band insulator is prepared, which has been 
achieved jl6j . It should be noted however, that current 
approaches to generating staggered magnetic fields with 
vector light shifts in the alkali's [TT] only work for the 
high-Z atoms Rb and Cs, for which there are no fermionic 
isotopes. On the other hand, spin dependent optical po- 
tentials are possible for alkaline earth atoms, including 
the fermionic isotopes [TT] . 

In principle, statements of the coherent quantum evo- 
lution may be more subtle than a question of the ground 
state. For example, the statement, that all energies get 
mapped to minus themselves upon mapping from % to 
— does not imply that expectation values of observ- 
ables evolve in exactly the same way as they would under 
the influence of —T-L. As we will now show, however, time 
reversal symmetry dictates that the evolution will indeed 
be the same in certain situations. 

We consider a general Hamiltonian assumed to 



be invariant under the application of some time reversal 
operator i?^, i.e., Rl'H{t)Rt = Hit). There is some free- 
dom in the choice of the time reversal operator. For our 
concrete example we take Rt to be the operator giving 
the complex conjugate in the basis {| t)j I 4-)} so that 
Rtia\ t) + b\ D) = a*\ t) + b*\ i). With this choice we 

have RjS^Rt = S'f , RlSfRt = Sf, and RjSfRt = -Sf 
so that the Hamiltonian (pi) is invariant under time rever- 
sal. Note that if the {| t): I i)} basis represents a physi- 
cal spin- 1/2 then the true time reversal is exp(i7rS'^)i?t. 
However, for our argument Rt need not be the physical 
time reversal operation. 

We now consider the time evolution of some operator 



{A{t))= (^o|T_i{e'^o'i*'«(*')} 
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where Tt (T^t) denotes (anti) time ordering. Since 
RtTte-'^fo dt'nt')j^l ^ y^g+./J dt'nt') 

it immediately follows that if we start in a state invariant 
under Rt {Rt\ipo) — IV'o)), then the time evolution of 
any invariant operator RjARt — A, will be the same as 
the time evolution under — ■H(t). For the preparation of 
AF both the initial state ([3| and the stagerred (Neel) 
magnetization are indeed invariant under Rt- 

The final Hamiltonian ([T]) for the AF states has a con- 
tinuum of (ungapped) excitations, which will always re- 
sult in some excitations for a finite preparation time. We 
analyze this imperfection in a three dimensional cubic op- 
tical lattices using a slightly modified version of standard 
spin wave theory. This approach is expected to give ac- 
curate result in the regime where the Neel magnetization 
is large, which is our main interest here. In the simu- 
lations we also include results far away from full mag- 
netization. These results are not quantitatively reliable 
but are indicative for the regime, where the preparation 
breaks down. 

To avoid dealing with two sublattices it is convenient 
to apply a mathematical transformation that rotates the 
spins of the B sublattice by tt around S^. In the ro- 
tated basis the Neel state is transformed to | tt ••• t)- 
Accordingly, the equations of motion become 



dt 



Ht)S+ (6) 



Here and below, " " " denotes the spin operators trans- 
formed by the unitary sub-lattice rotation. To find an 
approximate solution for the operators, we assume that 
in products of operators such as S^Sj we may replace 

Sj by its mean value. The mean spin is then determined 
self-consistently by requiring {Sp = 1/2 — {S~Sj~). 
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To describe the bulk properties of the system we as- 
sume periodic boundary conditions and switch to the mo- 
mentum representation 5'j^and = S^^ , where k is a 
reciprocal lattice vector. The equation of motion then 
only couples and SZ]^ and may be conveniently solved 
by a Bogoliubov transformation S^{t) ~ Ui^{t)S^{t = 
0) -t- Vk{t)SZ]^{t — 0). To allow for imperfections in the 
initial preparation of the atomic spin state we assume 
that each atom is prepared in the wrong spin state with 
a probability P. For small P we can represent this as a 
small contribution to all modes in the momentum repre- 
sentation and the mean value of the (Neel) magnetization 
is then given by 

i^') = l-p-^T.\-^\'- (7) 

k 

In Fig. [TJwe show the time evolution of the staggered 
magnetization obtained by numerically solving the equa- 
tions for Mk and Wk with the mean magnetization given 
by Eq. (j?]) on a lattice of size 45'^. For comparison we 
also show the magnetization in the instantaneous ground 
state found by diagonalizing the equation of motion for 
a fixed (-S*^), evaluating (S^) with Eq. ([t]), and iterat- 
ing the solution to convergence. In the figure we have 
used a magnetic field of the form h{t) — /iqC^"*, and we 
have excluded the contribution from the k = Goldstone 
mode since this mode only represent a slow rotation of 
the broken symmetry axis. Experimentally this mode 
may be effectively suppressed by not reducing the stag- 
gered magnetic field all the way to zero. 

As opposed to the ideal (classical) Neel state with 
(5^) = 0.5, the AF eigenstate of Eq. ^ has a slightly 
reduced value of (S^) due to quantum fluctuations as in- 
dicated by the dots in Fig. [l] The imperfection in the 
preparation is therefore characterized by the difference 
in (S^) between the ground and the prepared states. As 
demonstrated in Fig. [l] a) the adiabatic method can in- 
deed prepare states which are very close to the AF ground 
state. That the prepared state really is a (meta)stable 
state of the AF interaction may be verified experimen- 
tally by comparing the evolution with the evolution from 
the state with the opposite magnetization, which will not 
be stable [dashed curve in Fig. [l]a)]. In Fig. |2]we show 
the dependence of the final magnetization (S^) on the 
parameters of the protocol. As shown in the figure the 
procedure is applicable even for moderate field strengths 
ho > 5J and rather fast extinction rates a < lOJ. Fur- 
thermore the procedure is also robust against imperfec- 
tions in the preparation of the initial Neel state. We 
emphasize, however, that the states prepared in this way 
are not exact equilibrium states of the Hamiltonian ([T]) 
and following the adiabatic evolution the state may relax 
to a state with a lower magnetization. When the mag- 
netization is near the value in the ideal ground state, as 
in Fig. ([T]) a), we expect that this will not change the 
magnetization significantly. 

Ref. [7j investigated an adiabatic preparation scheme 




FIG. 1: Time evolution of the Neel magnetization when the 
staggered magnetic field is turned off. The full curves show 
the evolution for various imperfections in the initial prepara- 
tion (counting from above P = 0, 0.05, 0.1, 0.15, 0.2). For 
comparison the dots show the magnetization of the ground 
state of the instantaneous Hamiltonian. a) Adiabatic turn off 
with a = J and ho — 20J. b) Without staggered magnetic 
field ho = 0. The dashed curve in a) show the evolution start- 
ing from a state with the spin aligned oppositely (Sz) = —1/2 
for P = 0. 



using a slow increase of the tunneling r with respect to 
the interaction strength [/, but without an external stag- 
gered field. In the Mott state the spin interaction is ef- 
fectively described by Eq. ^ with J r^/C/. The 
unitary time evolution exp(— i J dt'H{t)) arising from a 
time dependent interaction strength in Eq. ^ is iden- 
tical to the evolution with an average value J given by 
Jt — J dt' J{t'). Adiabatic increase of the tunneling can 
therefore be mapped to our scheme, but with h{t) — 
throughout. The result of such a preparation scheme is 
shown in Fig. [l] b) , and this procedure always prepares 
states, which are far from the ground state. 

Our analysis so far assumed a truly bounded spectrum 
and therefore a state prepared within the highest energy 
manifold can in principle be infinitely long lived. In prac- 
tice, there are decay channels due to the existence of even 
higher energy states, which are not included in the effec- 
tive spin model ([T]). We should therefore ensure that the 
state docs not have time to decay during the preparation 
stage. 




12 12 
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FIG. 2: Neel magnetization {S^) at the end of the adiabatic 
evolution, a) As a function of the initial magnetic field ho 
and extinction rate a for a perfect initial state P = 0, and b) 
as function of extinction rate a and imperfection in the initial 
state P with a fixed initial field ho = lOJ. 
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One decay channel is provided by the soft compress- 
ible edges, which surround the Mott regions for inhomo- 
geneous trapped systems. Perturbations from the edge 
enter the bulk with velocity smaller than the spin wave 
velocity vs ~ J- To avoid ruining our state by pertur- 
bations coming from the edge, we need to be fast on a 
time scale given by todgo ~ N,n/vs, where Nm is the 
length of the Mott state. From Figs, l] a) and [2j we see 
that the preparation can be accomplished on a time scale 
^adiabat ^ 5/ J. If the Mott regions are sufficiently large 
so that Njn ^ 5, decay of the staggered magnetization in 
the bulk into edge modes is effectively suppressed during 
the preparation. 

Another relaxation mechanism owes to the existence 
of bulk particle hole excitations in the Mott state on 
which the effective spin Hamiltonian ([I]) is built. The 
high energy AF spin state can release energy into spin 
de-excitations while exciting a doubly occupied site to 
balance the energy cost. Since a single spin flip only car- 
ries away an energy of order J — At^/U, while the dou- 
ble occupation costs energy U, the decay can only occur 
in a rather high order process whereby a large number 
n ~ {U /t)"^ of spin de-excitations are created. A similar 
relaxation process was recently considered in analyzing 
the decay of "doublon" excitations in the Fermionic Hub- 
bard model IS]. The decay rate computed by Fermi's 
golden rule is small, and scales as F ^ t exp(— y4f7^/T^), 
where A depends only logarithmically on U /r and can be 
taken as a constant of order one for the experimentally 
relevant parameter regime U > lOr. 

A different imperfection comes from holes in the Mott 
state. Initially, nearest neighbor hopping of the holes 
will be completely suppressed by the strong linear con- 
finement in presence of the staggered field if r <C /iq- The 
only allowed hopping will be via a second order process 
and the initial state is thus a low (high) energy state of 
the AF (ferromagnetic) Hamiltonian even in the presence 
of holes. When decreasing the staggered field one adia- 
batically prepares low (high) energy states of the holes, 
which are not disruptive to the AF order. Thus, at least 
for low density of holes we do not expect them to disrupt 
the AF order. 

Before concluding we point out that the method con- 



sidered in this letter is not confined to spin systems. In 
particular there is an interesting application for interact- 
ing spinless bosons in an optical lattice. For usual boson 
hopping —JJ2{ij) ^i^j with J > 0, the hopping Hamilto- 
nian has a unique ground state in which the single par- 
ticle wave-function has a uniform phase. So, regardless 
of the lattice geometry, a macroscopic number of bosons 
would simply condense in this state. 

The Hamiltonian with J < is more interesting. The 
single particle ground state seeks an optimal phase dif- 
ference of TT across all links, which is not always possible. 
Certain lattice geometries, are then frustrated because 
they do not support a unique ground state phase config- 
uration. In the Kagome lattice, for example, the effect 
is particularly dramatic, resulting in a degenerate band 
of single particle ground states [H]. Because the kinetic 
energy is completely quenched, the many-body ground 
state is crucially determined by interactions, regardless 
of how weak they are. How exactly the quantum frus- 
tration is relieved is not known, but the process is likely 
to give rise to an exotic phase. The general scheme we 
propose provides a natural route to realize such states as 
the highest energy states of a lattice boson Hamiltonian 
with usual hopping and attractive interactions between 
bosons. Then, —H will implement positive hopping and 
repulsive interactions between particles. 

In conclusion we have shown that adiabatic prepara- 
tion allows the study of interesting many body dynamics 
as the highest energy state of Hamiltonians whose ground 
states are less interesting. As a particular example we 
have analyzed in detail the preparation of AF states us- 
ing this approach. The method is argued to be robust 
to defects in the initial configuration and allows for fast 
preparation. We argued that the same approach can be 
used to prepare interesting states of lattice bosons with 
a geometrically frustrated hopping term. 
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